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Abstract 

Energetic positrons produced in annihilation or decay of dark matter particles in the 
Milky Way can serve as an important indirect signature of dark matter. Computing 
the positron flux expected in a given dark matter model involves solving transport 
equations, which account for interaction of positrons with matter and galactic mag- 
netic fields. Existing calculations solve the equations inside the diffusion zone, where 
galactic magnetic fields confine positrons, and assume vanishing positron density on 
the boundaries of this zone. However, in many models, a substantial fraction of the 
dark matter halo lies outside the diffusion zone. Positrons produced there can then 
enter the diffusion zone and get trapped, potentially reaching the Earth and increasing 
the expected flux. We calculate this enhancement for a variety of models. We also 
evaluate the expected enhancement of the flux of energetic photons produced by the 
inverse Compton scattering of the extra positrons on starlight and cosmic microwave 
background. We find maximal flux enhancements of order 20% in both cases. 



1 Introduction 



Multiple observations, ranging from rotation curves of galaxies to the structure of anisotropies 
in the cosmic microwave background, indicate the presence of substantial amount of dark 
matter in the universe. While the microscopic nature of dark matter has not been probed so 
far, one popular hypothesis states that dark matter consists of weakly interacting massive 
particles (WIMPs). The relic density of WlMPs, predicted within the simplest thermal- relic 
cosmological scenario, turns out to be consistent with the observed dark matter abundance. 
This coincidence, and the fact that WIMPs are in fact predicted by many popular extensions 
of the standard model of particle physics, make the WIMP hypothesis quite attractive on 
theoretical grounds. 

If dark matter in the Milky Way halo indeed consists of WIMPs, occasional pair-annihilation 
of these particles should produce spectacular high-energy cosmic rays, giving a potential 
indirect signature for dark matter. In particular, energetic positrons produced in WIMP an- 
nihilation events provide a promising way to look for dark matter, which has been exploited 
by a number of experiments. Two satellite-borne experiments, HEAT [1] and PAMELA [2], 
provided measurements of the positron flux in the 1 — 100 GeV range. In addition, several 
recent experiments - most notably ATIC [3], FERMI [4], and HESS [5] - have measured the 
sum of electron and positron fluxes, extending to TeV energies. All of these experiments re- 
port fluxes and spectra inconsistent with "canonical" background models. The excess fluxes 
can be interpreted as coming from dark matter annihilation; while the minimal WIMP mod- 
els tend to predict fluxes too low to be consistent with experiments, simple extensions of the 
WIMP paradigm, e.g. incorporating the Sommerfeld mechanism to enhance the annihilation 
rates at low velocities [6] , can flt the data. A decaying dark matter particle with lifetime of 
order 10^^ sec is another candidate [7]. It should be noted that many uncertainties remain 
in the evaluation of positron fluxes from conventional astrophysical sources, and viable ex- 
planations of the observed excesses in terms of conventional astrophysical sources have been 
proposed [8, 9]. Still, the recent rapid experimental progress in this area has highlighted the 
need for accurate predictions of positron fluxes, both from conventional and unconventional 
{e.g. dark matter) sources. 

Our focus in this paper is on the calculation of positron fluxes from dark matter annihila- 
tion. The calculation proceeds in two steps. First, the spectrum of positrons emerging from a 
WIMP pair-annihilation event is calculated. (This spectrum depends on the particle physics 
model responsible for WIMPs: for example, WIMPs can annihilate directly into e+e~ pairs, 
or into W'^W~ with subsequent decay — >■ e+i^e, etc.) Second, the interactions of the 
positrons with galactic magnetic flelds, starlight and CMB photons, synchrotron radiation, 
and other effects occurring on the way from the production point to the detector must be 
included. The propagation of positrons through the galactic medium is governed by the 
transport equations, generally a complicated system of coupled differential equations involv- 
ing densities of positrons, photons and other cosmic ray species. Comprehensive, detailed 
description of galactic propagation requires numerical techniques, and extensive packages 
such as GALPROP [10] have been developed to tackle the problem. However, a reasonable flrst 
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approximation to positron propagation can be obtained by treating the positron density in 
isolation, and modeling positron interactions with the medium by simple linear diffusion and 
energy-loss terms (see e.g. Refs. [11, 12]). This results in the "diffusion- loss" equation: 



where ipijKjEjt) = dn^.+ /dE is the positron density per unit volume per unit energy; K is 
the diffusion coefficient describing the interaction of the positron with the stochastic galactic 
magnetic field; b is the rate at which positrons lose energy due to synchrotron emission and 
inverse Compton scattering (ICS); and q is the source term, in our case due to dark matter 
annihilation or decay. It is reasonable to assume that the present positron density reflects 
a time-independent steady state, and solve Eq. (1) with dip/dt = 0. To do this, boundary 
conditions in space need to be imposed. Galactic magnetic fields are confined to a cylindrical 
"diffusion region" or "diffusion zone", with radius R of order 20 kpc and half-thickness L 
taken to be between 1 and 15 kpc. If a positron is injected inside the diffusion region, it slowly 
random- walks through the region, taking of order 10^ years (for typical parameters) to reach 
the edge. Once the edge is reached, the positron escapes into intergalactic medium, traveling 
essentially with the speed of light. Thus, in steady state, the positron density outside the 
diffusion region is expected to be strongly suppressed, and flux calculations assume = at 
the diffusion region boundaries, z — ±L and r — R. The transport equation is then solved 
within this cylindrical region. Essentially all existing calculations of positron flux from dark 
matter annihilation, analytic or numerical, make this assumption.^ The key observation of 
this paper is that this choice of boundary conditions results in an under-estimate of the 
positron fluxes from dark matter annihilation. 

The source of the flux enhancement is illustrated in Fig. 1. A typical dark matter 
halo is spherically symmetric and extends beyond the diffusion region, in particular in the 
vertical direction: for example, for an isothermal dark matter profile in the M2 propagation 
model (which uses L = 1 kpc), the diffusion zone contains only 10% of the dark matter 
mass of the full halo. An order-one fraction of the positrons produced by dark matter 
annihilations outside the diffusion region will enter the diffusion region, and get "stuck" 
there. These positrons will contribute to the steady-state density (and fiux) inside the 
diffusion region. This contribution is missing from any calculation that only considers the 
sources at \z\ < L. It cannot be incorporated by simply extending L to a larger value since, 
for a given propagation model, L is a fixed physical parameter that represents how far the 
galactic magnetic fields extend in space. What is needed is a formalism that keeps galactic 
magnetic fields confined within the region defined by L as required by the propagation 
model being used, yet incorporates contributions from dark matter annihilation beyond this 
region. In this paper we present such a formalism, via a simple extension of the Bessel- 
transform approach of Ref. [12]. We then analyze its impact quantitatively, in a number of 

^Some of the uncertainties in dark matter indirect signals associated with modeling the diffusion zone 
boundary have been recently studied in Ref. [13] , using the model with a position-dependent (exponential in 
z) diffusion coefficient proposed in [14]. Also, the positron flux from galactic subhalos has been considered 
in Ref. [15]. 
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Figure 1: Left - The dark matter halo extends significantly beyond the diffusion zone, but 
only sources inside the zone are considered in the conventional formalism (Section 2.1). 
Right - The extended formalism (Section 2.2) considers sources in the free propagation zone 
in addition to the diffusion zone. 

illustrative models. We also analyze the contribution of the extra positrons to the expected 
flux of energetic photons from the inverse Compton scattering (ICS) process in the galaxy, 
suggested in Ref. [16] as a robust signature of strongly annihilating dark matter. 

2 Formalism 

In this section, we will first review the conventional Bessel-transform approach [12] to solving 
Eq. (1), and then present our extension of this formalism. 

2.1 Conventional Formalism 

In the conventional formalism, one assumes that the diffusion coefficient and the energy loss 
term are position-independent within the diffusion zone, but depend on energy. To describe 
the energy dependence, we adopt the simple power-law parametrization of Ref. [12]: 

KiE) = Koe', h{E) = ^e\ (2) 

where e = E/Eq, Eq = 1 GeV. The quantities Kq, 6, and te, along with R and L, define 
the "galactic propagation model". Standard choices are te = 10^^ sec and R = 20 kpc. 
We will use several combinations of the other parameters, summarized in Table 2. These 
combinations have been found compatible with observed cosmic ray properties, such as the 
B/C ratio [17]. 

Assuming that the positron source term is cylindrically symmetric, and imposing the 
boundary condition ip = at r = R and z = ±L, the solution to Eq. (1) can be expressed 
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as a Bessel-Fourier series: 
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where Jj denotes the i-th order Bessel function of the first kind, and a^'s are the zeros of Jq. 
Taking the Bessel and Fourier transforms of Eq. (1) at steady state and changing variables 
according to 
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where the coefficients Qi^n are the Bessel-Fourier transforms of the source term: 
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In dark matter applications, the positron source energy spectrum is position- independent: 
q{r, z, e) = R{r, z)f{e), where /(e) = dNe+ / dEe+ is the energy distribution of positrons from 
a single WIMP annihilation or decay. Specifically, for annihilating dark matter, R{r, z) = 

ri{av) 



"'~x 



, where rj is equal to 1/2 for Majorana WIMPs and 1/4 for Dirac WIMPs, p 
is dark matter energy density, is the WIMP mass, and {av) is the thermally averaged 
cross section for WIMPs to annihilate into a final state containing positrons. For decaying 
dark matter, R{r,z) = Vp{r,z)/m^, where F is the partial decay width into a final state 
containing positrons. For such sources, the positron density can be expressed as 
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where Cmax is the maximum energy at which positrons are produced (typically equal to m-^) , 
and 
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is the halo function. Here Ri^n are coefficients of tfie Bessel- Fourier expansion of the function 
R{r, z) . Note that the halo function only depends on e and es through the combination 



where Ad is the diffusion length. 
2.2 Extended Formalism 

To include contributions from dark matter annihilations or decays outside the diffusion zone, 
we would like to solve Eq. (1) in a larger cylinder extending out to \z\ = D, with D > L such 
that all (or essentially all) of the dark matter is contained within this cylinder. We impose 
the boundary condition if) — Q &X, z — ±£), since there is no influx of positrons through the 
boundaries. Inside the cylinder, there are two zones: diffusion zone., extending out to |2;| = L, 
and free-propagation zone., with L <\z\ < D. We expect that the positron density is small 
in the free-propagation zone, since the positrons propagate through this zone very quickly; 
however, there is an influx of positrons from the free-propagation zone into the diffusion zone, 
which increases the steady state flux there. For mathematical convenience, we model the 
propagation through the free-propagation zone by the same diffusion-loss equation, Eq. (1), 
but with a different diffusion coefficient = Kie^ . Ideally we would like to take the limit 
Ki — 7- oo, corresponding to very fast diffusion. In practice, numerical convergence issues put 
an upper limit on the ratio of Ki/Kq; however, we are able to calculate with values of Ki for 
which the results are essentially independent of that parameter, suggesting that the limiting 
behavior has been reached. In principle, wc would also like to set the energy loss term to 
zero in the free-propagation zone, 6 — )> 0. However, since the positrons spend negligible time 
in this zone, in practice this term does not affect the results, and for simplicity we use the 
same b in the diffusion and free-propagation zones. 

With a position-dependent diffusion coefficient, the steady state transport equation is 

-VK -V^-KAi;- A[ft(x, E)^;] = g(x, E) . (12) 

Since both K and q have cylindrical symmetry, we again expand ijj and q in Bessel-Fourier 
series as in Eqs. (3) and (6), but with L ^ D. Substituting these expansions into Eq. (12) 
yields (for each i): 

.(.e)I:A.(e)((^)^(i)>.(^) 
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In our setup, the diffusion coefficient lias tfie form 



K{z,e) = (Ko + k{z))e\ (14) 

wliere X — > in tlie diffusion zone and K ^ Ki — Kq Ki Ki in the free-propagation zone. 
Multiplying both sides of Eq. (13) by sin ( ^^^2d^^ ^ (where m is an integer) and integrating 
over z e [— -D, D] gives 

/ nvr \ dK (mx{z + D)\ (m-n{z + D)\ , 

( ) j ^i,m + 'dt^^'^ ~ Qi,m , (15) 

where we performed a change of variables as in Eq. (4). The last line is just what one would 
obtain in the conventional formalism with L ^ D; however, there are now additional terms 
that mix different Fourier components of the positron density. It is useful to put this system 
of equations in matrix form: 

dP- 

y + A,-P, = Q,, (16) 

where Pj and Qj are vectors containing the Pj^„ and Qi^n, and Aj are matrices whose elements 
can be read from Eq. (15). The matrices Aj are f-independent. The solution is given by 




Pi{t) = / dts exp[-(i - QA^] Q, . (17) 

The positron density from dark matter decay or annihilation has the same form as in the 
conventional formalism, Eq. (9), with the halo function (10) replaced by 

I{r,z,e,es) = EE (^) sin ( ^^^^ + (exp[-(t - t,)A.]R.)^ . (18) 

As before, the halo function only depends on e and through t — ts', however, in this case, 
the definition of diffusion length A^; is ambiguous, since it involves the diffusion coefficient 
which is now ^-dependent. Below, we will always use the value of the diffusion coefficient 
inside the diffusion zone to define Ad, in complete analogy with Eq. (11). 

To proceed with the analysis, we need to specify precisely how the diffusion coefficient 
depends on z. The simplest choice is to model it as a step function with the conventional 
value Kq in the diffusion zone and a much larger value Ki in the free-propagation zone. Such 
a model, however, is ill-suited for the Bessel- Fourier approach, since an infinitely sharp jump 
in K requires a very large number of terms in the expansion to achieve convergence. To 
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Model 


a 




7 


rs(kpc) 


Cored isothermal 


2 


2 





5 


NFW 


1 


3 


1 


20 


Moore 


1.5 


3 


1.3 


30 



Table 1: Dark matter density distribution profiles. 



avoid this problem, we smooth out the step function across a finite interval L < z < L -\- 
with d <^ L,D. Specifically, we assume 



K{z) = 



\{K^ + Ko) -\{K^- K,) [3 cos (^tt) - cos^ (^tt) 



if 1^1 < L; 

ii L < \z\ < L + d; 



iiL + d< 



<D. 

(19) 

This function has continuous first and second derivatives at z = L and z = L + d. Another 
important advantage of this form is that all integrals in Eq. (15) can be evaluated analytically, 
leading to significant speedup of numerical calculations. It should be kept in mind that while 
in reality the galactic magnetic fields probably do drop off smoothly over some finite distance 
at the boundary of the diffusion zone, the particular choices of d and the analytic form of 
K{z) that we make are not physically motivated. We will show that the results of our 
analysis are approximately independent of these choices. 



3 Dark Matter and Galactic Propagation Models 

The dark matter density distribution in the Milky Way halo is modeled with the generic 
profile 

where = 8.5 kpc is the distance from the solar system to the galactic centre, and Pq = 
0.3 GcV cm~^ is the local dark matter density in the solar neighborhood. Wc use three 
profiles: isothermal, Moore, and Navarro, Frenk and White (NFW) (see Table 1). Since the 
numerically derived NFW and Moore profiles diverge at the center of the galaxy, the profile 
inside r< tq is replaced by the smoother profile, as in Ref. [12]: 

p(r) = po(l + fli sinh(7ra;) + 02 sinh(27ra;))^^^ , (21) 

where X = r/ro, po = p^ro), ai = 02 + 27, and 02 = 87(7r^ — 9 + 67))/(9(3 — 27)). This 
rcnormalized profile has a continuous first derivative at ro and preserves the total number of 
annihilations within the core. 
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Model 


S 


Ko (kpcVMyr) 


L (kpc) 


MED 


0.70 


0.0112 


4 


M2 


0.55 


0.00595 


1 



Table 2: Galactic propagation models. 



The positron injection spectrum f{E) depends on the microscopic model of dark matter, 
in particular the WIMP mass and its annihilation channels or decay pattern. Motivated 
by PAMELA and FERMI data, we consider a rather heavy WIMP, = 3 TeV (see, 
e.g., Ref. [18]). Since current data favor leptophilic models {i.e. those where dark matter 
annihilations or decays result mostly in all-leptonic final states), we focus on the following 
four annihilation channels: 

1. XX ~^ e+e". This channel is not favored by PAMELA and FERMI data [18], but 
produces a simple monochromatic injection spectrum and hence serves as a useful 
limit. 

2- XX ~^ A*^/^") with the muons then decaying to give positrons via the familiar decay 

3- XX — >■ 00 — >■ 4e, where (p is some intermediate particle of mass (for example, a 
gauge boson of an extra gauge symmetry that dark matter is charged under, as in 
Refs. [6, 20] and others). 

4- XX — > 00 — > 4/i. The motivation is the same as in 3, but the positron injection 
spectrum is softer in this case. 

It is straightforward to obtain f{E) for each of these models. We consider the same final 
states for the case of decaying dark matter, and in fact the only difference is that the energy 
of the primaries (e in process 1, in process 2, particles in processes 3 and 4) is halved 
compared to the annihilating DM scenario. 

For positron propagation in the galaxy we use the M2 and MED propagation models [12] 
(see Table 2), which are compatible with cosmic ray data. We expect the enhancement of 
the positron fiux at the solar position to be most significant in models with small L, since 
this both maximizes the amount of dark matter outside the diffusion zone and minimizes the 
energy losses of the positrons coming from outside the diffusion zone. Because of this, we 
do not study other well-known galactic propagation models such as Ml or MAX, in which 
L > 10 kpc and no significant fiux enhancement is expected. 

To complete the description of galactic propagation in the extended formahsm of Sec- 
tion 2.2, we need to specify the parameters D, Ki and d; these, and some other numerical 
issues, are discussed in the Appendix. 
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Figure 2: Examples of conventional (dashed) and extended (solid) halo functions. 

4 Results: Positron Fluxes 

We analyze the two galactic propagation models, three galactic halo profiles, and four dark 
matter annihilation scenarios specified in Section 3. In each model, we computed the halo 
function (18) using the extended formalism of Section 2.2, which was then used to compute 
positron density at the solar position (r© = 8.5 kpc, Zq = 0) as a function of energy. For 
comparison, we also computed the halo function and positron density within the conventional 
formalism (Sec. 2.1), which neglects the contribution due to dark matter annihilations outside 
of the diffusion zone. The positron flux measured by a Solar System based experiment is 
given by 

^AE) = ^^{r^,Ze,E) (22) 

where /3e+ is the velocity of a positron of energy E. 

Figure 2 shows the correction to the halo function for a few cases. The corrections are 
generally small, even for the M2 model with L — 1 kpc. As expected, the correction is larger 
for flatter profiles and in regions closer to the diffusion zone boundary. 

In Fig. 3, we plot the ratios of the full positron flux and the flux computed within the 
conventional formalism for various scenarios. Overall, the enhancements we observe are not 
large, typically in a few-% range, and are thus smaller than other astrophysical uncertainties 
at this point. The plots in Fig. 3 demonstrate many expected features. First, the enhance- 
ment is largest for annihilation into e^e~, which has the most energetic input spectrum, and 
progressively decreases for less energetic input spectra. Second, the plots are consistent with 
the notion that for M2 propagation the spectrum at Earth is influenced largely by sources 
within a few kpc and, in particular, is insensitive to the large cusp at the center of the 
galaxy. In particular, although the Moore profile is more cusped at the galactic center than 
the NEW profile, in the solar neighborhood the dark matter density drops off with z faster 
in the NEW profile than in the Moore profile. Consistent with this, the fiux enhancement 
is smaller in the NEW profile. Third, the enhancement is larger at lower positron energies, 
consistent with the fact that positrons lose energy as they propagate through the intergalac- 
tic medium and hence energetic positrons entering from the halo outside the diffusion zone 
will arrive at the Earth at lower energies. Fourth, the enhancement is almost negligible 
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Figure 3: Positron flux enhancement for annihilating dark matter, for = 3 TeV. 



for MED propagation, since there is negligible amount of dark matter outside the diffusion 
zone compared to inside it in this case, and the diffusion zone boundary is also farther away 
from the Earth. All these features confirm that the extended formalism we suggest, and our 
numerical approximations, are physically sensible. 

4.1 Decaying Dark Matter 

Dark matter decaying with a lifetime of about 10^ times longer than the age of the universe 
can also explain the PAMELA and FERMI excesses [7]. While dark matter annihilation rates 
are proportional to p^, the rate for decaying dark matter is proportional to instead, which 
implies that the relative contribution from the halo exterior to the diffusion zone should be 
greater. 

We computed the enhancement for decaying dark matter for the isothermal profile since, 
for a given galactic propagation model, the enhancements are the largest for this profile. 
The results are plotted in Fig. 4 and should be contrasted with the bottom row of Fig. 3. 
For MED propagation the enhancement remains negligible; for M2 propagation the increase 
over the annihilating dark matter scenario is very small because the region of influence only 
extends to a few kpc, where differences between (p^) and {p^) are small. 
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Figure 4: Positron flux enhancement for decaying dark matter, for m^ = 6 TeV. 

5 Gamma Rays from Inverse Compton Scattering 

Tlie positron flux corrections calculated in the previous section are at the solar position, 
which lies at 2; = and hence only gets minimal corrections. The more significant corrections 
occur close to the diffusion zone boundary, but these cannot be observed directly. On 
the other hand, gamma ray flux from inverse Compton scattering (ICS) off positrons is 
sensitive to positron density throughout the galaxy, since a photon scattering off an energetic 
positron towards the Earth anywhere in the galaxy will travel unperturbed through the 
interstellar medium and can be detected. Hence large corrections to positron density close 
to the boundary can significantly affect the ICS spectrum measured at Earth. 

A semi-analytic calculation of the ICS energy spectrum from dark matter annihilation is 
presented in [18] and [19]. The ICS flux is expressed (see Eq.(lO) of [18]) as 

c[$y^^ 9rQ{av) / Pq\ 

Here we have combined the two separate dimensionless parameters Jjc and Gjc defined in 
Eq. (11) of Ref. [18] into a single parameter 

Huc = mt fdnf ^U{E,) [ rdE'f{E')m,E',r). (24) 

J Jl.o.s Vq Utot J E^ J El R[Ee) JEe 

This is necessitated by the explicit position dependence of the halo function / in our treat- 
ment, which requires it to be integrated over both position and energies. E^ and Ey denote 
photon energies before and after scattering. The sum over i accounts for the three compo- 
nents of galactic light that can scatter off energetic positrons: CMB, starlight, and starlight 
rescattered by dust. The energy density profiles Ui, energy spectra f^i, mean energies {E^i), 
the inverse Compton factor fjc in the scattered photon spectrum, and relativistic correc- 
tion R{Ee) are as presented in Ref. [18]. Eq. (24) shows that if the line of sight ends at 
the diffusion zone boundary, the halo function I at that position, which receives significant 
corrections from the extended halo, directly enters the calculation for the ICS flux. 
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Figure 5: Left column: Halo function from conventional (dashed) and extended (solid) 
formalisms as a function of position along line of sight for Xd = 0.5 kpc. Right column: ICS 
gamma ray flux enhancement for annihilating dark matter, = 3 TeV. 



As an example, we consider the line of sight from the Earth (r = 8.5 kpc, z = 0) 
to the edge of the diffusion zone closest to the galactic center {r = 0, z = 1 kpc in the 
M2 model). In this case, the largest contribution to the ICS flux enhancement comes from 
regions close to the diffusion zone boundary: the photon density is the greatest there because 
of the proximity to the galactic center, and the positron density is zero or close to zero in 
the conventional formalism because of the boundary conditions but can be significant in 
the extended formalism (see Fig. 5, left column). The second column in Fig. 5 shows the 
enhancement in the ICS flux, calculated using Eq. (23) and Eq. (24). As expected, the 
enhancement in the ICS gamma ray flux is more significant than in positron flux. While 
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there is essentially no correction to the positron flux around the TeV scale, the ICS flux at 
this energy receives corrections of about 20%. Note that the enhancement plots - especially 
in the e^e~ channel - show three distinct "bumps". These correspond to the three different 
galactic gamma-ray components: CMB, starlight, and starlight rescattered off dust. 

6 Conclusions 

We summarize our flndings as follows: 

• Including contributions from the extended dark matter halo to positron density in 
the galactic diffusion zone coming from dark matter annihilation or decay can result 
in corrections of up to 17% in the observed positron flux and up to 24% in the ICS 
gamma ray flux. 

• For positron flux, the enhancement is progressively lower at higher energies since 
positrons coming in from the halo lose energy in the diffusion zone and arrive at 
detectors with lowered energies. While we see enhancements of up to 17% at 10 GeV, 
enhancements are <10% for all considered cases for energies >100 GeV. 

• For ICS gamma ray ffux, the bulk of the correction comes from photons scattering 
off energetic positrons close to the diffusion zone boundary, hence the enhancement is 
maintained even at the highest photon energies. 

• The enhancements are most signiflcant for the M2 propagation model, where the diffu- 
sion zone is only 2 kpc thick and hence the bulk of the dark matter halo lies outside this 
zone. For MED model (and presumably other models with larger L), the corrections 
are negligible. 

• At present, experimental uncertainties on flux measurements as well as astrophysical 

uncertainties in the positron and ICS fluxes - these come from numerous sources, such 
as uncertainties from dark matter profiles, propagation models, energy density and 
spectra of photons in the galaxy, and the simplifications made to the transport equation 
to describe positron propagation - remain significantly greater than the additional 
contribution from the dark matter halo beyond the diffusion zone. It does not need to 
be included in fits to data at this stage, but should be considered when accuracy to 
better than 25% is needed. 

The extended formalism can also be augmented in a straightforward manner to include 
contributions from other important sources, such as dwarf galaxies, that are impossible to 
incorporate in the conventional formalism. 
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A Comments on Numerical Issues 



This appendix discusses some details of our numerical calculations. 

Choice of parameters — Parameter choices were dictated by the need to carry out all 
computations in a reasonable amount of time. For each dark matter density profile, D was 
chosen such that at (r = 0, z = D) is 5% of the corresponding value at the diffusion zone 
boundary (r = 0, 2; = L). Ki was typically chosen to be 2000 times ifo, while we typically 
used d = 0.2 kpc and d = 0.5 kpc for M2 and MED propagation respectively. Although 
R = 20 kpc is the standard radius of the diffusion zone cylinder, we used R = 11 kpc for runs 
with NFW and Moore profiles to save computation time. This choice is justified because 
(i) there is negligible amount of dark matter beyond iZ = 11 kpc for these profiles, and (ii) 
all our results for these profiles use the M2 propagation model, where the diffusion zone 
cylinder height is so small that positron abundance is determined primarily by losses in the 
vertical (z) direction, and the diffusion region forgone by using the smaller radius does not 
significantly affect the results. 

The results should not be sensitive to the choices of Ki or d. To verify this, we varied Ki 
and d and checked how this affected our results. This is shown in Fig. 6 for the isothermal 
profile, M2 propagation. The first column shows the effects on the positron fiux at Earth. 
The top plot shows how the halo function changes at the solar position: the dashed and solid 
curves correspond to halo functions in the conventional and extended formalisms, while the 
two curves in between correspond to d decreased by 40% (dot dashed curve) and Ki increased 
by 45% (dotted curve) respectively. The next two plots show the corrections to the positron 
flux from these two parameter variations; the correction factor on the y-axis is calculated 
as the ratio of the flux after parameter variation to the flux before parameter variation. 
Likewise, the second column shows the corresponding effects on the ICS spectrum. The top 
plot shows how the halo function changes (along the line of sight used in Section 5) near the 
diffusion zone boundary, where the effect should be the most important, for A/j = 0.08 kpc 
(right). The next two plots show the corrections to the ICS flux from the two aforementioned 
variations. The sensitivity to these variations is <3%, while the effects we are studying give 
enhancements of 10-20% (Fig. 3 and Fig. 5). Therefore sensitivity to unphysical parameter 
choices is small compared to the physical correction we are studying, and does not affect out 
conclusions. 

To get results to converge, 2500 terms were needed in the Fourier series expansion. The 
Bessel series expansion required 60 terms for the isothermal proflle and 625 terms for the 
NFW and Moore profiles. 

Evaluating the halo function — To evaluate the exponential in Eq. (18), we diagonalized 
the Aj matrices. In a basis Vj in which Dj = V^^AjVj is diagonal, Eq. (18) becomes 

I{r,z,e,es) = E E -^o (^) sin (^^^^i|±^ j (v,exp[-(t - ^5)0,] Vr^R,)^ , (25) 

and the exponential is straightforward to calculate. This method was found to be more 
efficient than Taylor-expanding the exponential. The diagonalization procedure for large Aj 
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Figure 6: Halo function, positron flux, and ICS flux sensitivity to changes in d and Ki. 



can be sped up by noting that its entries are largest along the diagonal and get progressively 
smaller away from the diagonal, so we can approximate Aj to be n x n block diagonal, 
hence splitting the diagonalization task into n smaller, independent tasks. The 2500 x 2500 
Aj matrices were approximated to be block diagonal, containing five 500 x 500 blocks; to 
minimize the error from this approximation on the edge terms of the blocks, each block was 
diagonalized by first diagonalizing a larger 600 x 600 block, and throwing away 50 terms 
from each end. The results from this procedure were found to be in agreement with results 
obtained by diagonalizing the full Aj matrix for a few cases that were checked. 

Boundary position correction — The form of the diffusion coefficient in Eq. (19) has many 
benefits, but it also has a flaw: we ideally want K{z) to quickly deviate from Kq immediately 
outside the boundary z = L, but the function in Eq. (19) rises slowly at flrst before becoming 
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steeper because of the way trigonometric functions are shaped. Thus, over a small distance e 
just outside the boundary, K{z) ~ Kq, and this region L < z < L + e traps positrons just like 
the diffusion zone. For a fairer comparison, it is therefore necessary to shift the diffusion zone 
boundary inwards by a distance e (or move the boundary for the conventional calculation 
outwards by distance e). We choose e such that K{L + e) = 1.25Ko; this corresponds to 
about 5% of the thickness d. The motivation for this choice is that solutions obtained by 
setting dark matter sources to zero outside L + e and using the extended formalism with 
the diffusion zone boundary at z = L approximately match the solutions obtained from the 
conventional formalism with the diffusion zone boundary at z = L + e, which suggests that 
the enhancements must then come from contributions from dark matter sources in the free 
propagation zone. 

Consistency check — As an overall consistency check, we also reconstructed the positron 
density (Eq. (3)) in our computation and verified that it satisfies the diffusion-loss equa- 
tion (1) throughout the diffusion zone for 6 = 0. This was done for all results discussed in 
this paper. 
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